Transport on exploding percolation clusters 
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We propose a simple generalization of the explosive percolation process [Achlioptas et at, Science 
323, 1453 (2009)], and investigate its structural and transport properties. In this model, at each 
step, a set of q unoccupied bonds is randomly chosen. Each of these bonds is then associated with a 
weight given by the product of the cluster sizes that they would potentially connect, and only that 
bond among the q-set which has the smallest weight becomes occupied. Our results indicate that, 
at criticality, all finite-size scaling exponents for the spanning cluster, the conducting backbone, 
the cutting bonds, and the global conductance of the system, change continuously and significantly 
with q. Surprisingly, we also observe that systems with intermediate values of q display the worst 
conductive performance. This is explained by the strong inhibition of loops in the spanning cluster, 
resulting in a substantially smaller associated conducting backbone. 
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The quest for a percolation paradigm that displays a 
first-order transition has been the focus of intensive re- 
search in Statistical Physics 0. Very recently, an ex- 
tension of the traditional percolation model has been 
proposed that displays a sudden transition of the order 
parameter as a function of the bond occupation proba- 
bility 0. In this so-called "explosive percolation" (EP) 
process, bonds are selected to be occupied in accordance 
with a product rule that favors the growth of smaller clus- 
ters over large ones. Interestingly, although this model 
presents an abrupt transition when applied to different 
network topologies 0-11]) its critical phase reveals signa- 
tures of a typical continuous transition, as for example, 
a power- law distribution of cluster sizes 0. This con- 
troversial nature of the EP transition has been recently 
addressed in Ref. 0, where it has been argued through 
analytical arguments that the transition may actually be 
continuous for the case of random graphs. As a potential 
application, the EP model has been linked to the growth 
dynamics of Protein Homology Networks 0. Abrupt 
transitions have also been recently observed for other 
percolation-like models 0, . For instance, it has been 
shown that, at least in the limit of infinite dimensionality, 
the absence of loops represents an important ingredient 
to obtain an explosive percolation process [loj . Subse- 
quently, a study on regular lattices in which the growth 
of the largest cluster is systematically controlled during 
the percolation process, provided unambiguous evidence 
for a first-order phase transition 

The problem of transport in disordered media has 
been extensively studied under the traditional paradigm 
of percolation theory [l2|, [r|. Metal-insulator transi- 
tions [ljj and flow through porous media [l5[ are just 
two among several physical phenomena that have been 
successfully modeled under the framework of percolation 
as a characteristic structural second-order phase transi- 
tion. At the critical point, percolation systems typically 
display a conductance a that decays with the lattice size 
as a ~ with an exponent [i that depends exclu- 

sively on the lattice dimensionality [IH, [l3j]. Here we 



investigate the structural and transport properties of EP 
networks built on a square lattice. In our simulations 
a generalization of the product rule proposed in Ref. 
is adopted, where at each step a bond is selected from a 
random set of q unoccupied bonds. Strikingly, our results 
suggest that the scaling exponents associated with struc- 
tural and transport quantities change continuously with 
the parameter q. We observe that the abrupt behavior 
of the percolation transition is stressed as q increases, 
reaching a limiting regime for large q values, similar to 
the model where the selected connection is the one with 
the smallest possible weight in the entire network [Tfj| . 
Moreover, the global conductivity of the system exhibits 
a very slow decrease with system size for large values of q, 
with an exponent that approaches zero, while at interme- 
diate values the exponent is maximum, i.e., it decreases 
faster. 

The present model is implemented on a square lattice 
of size L. Initially all the bonds of the lattice are removed 
and then occupied in sequential steps. At each step, a set 
of q unoccupied bonds is randomly chosen. To each of 
these bonds a weight is assigned that is proportional to 
the product of the size (number of sites) of the clusters 
it would potentially join [2j. In the case a bond connects 
two sites in the same cluster, the weight is equal to the 
square of the cluster size. From the random set of bonds 
chosen, only the bond which has the smallest weight is 
occupied. The remaining bonds stay unoccupied, but 
can be selected again in later steps. For q = 1 we re- 
cover traditional percolation, while q — 2 corresponds 
to the Achlioptas model [2j. We stop the process when 
a spanning cluster is formed that connects the top and 
bottom of the lattice. At this point we apply the burning 
method [l7j to compute the mass of the spanning clus- 
ter M c i us , the mass of the conducting backbone M& oc fc, 
and the number M cut of cutting bonds (bonds that would 
disconnect the spanning cluster and stop conduction, if 
removed). The global conductance a of the system is 
calculated by solving the set of linear equations given by 
KirchhofFs conservation law on each site. Each of these 
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FIG. 1: (color online) Typical clusters obtained from the EP percolation process with g-links on a square lattice of size L = 512. 
In this model, at each step a set of q unoccupied bonds is randomly chosen. Each of these bonds has a weight given by the 
product of the cluster sizes that they would potentially connect. Only the bond with minimum weight in the set becomes 
occupied. The process continues in sequential steps, and is halted when a cluster forms connecting the top and the bottom 
of the lattice. Here we show clusters for q — 1 (traditional bond percolation model), q = 2 (Achlioptas model fjj), and 
q — 10. The spanning cluster appears in gray, the conducting backbone in black, and the cutting bonds in red. For q = 10 
the spanning cluster occupies almost the entire lattice and the conducting backbone constitutes a large fraction of the whole 
system. Interestingly, the number of cutting bonds is significantly larger in the case q — 2. 



quantities are averaged over at least 2000 realizations for 
each value of L. 

Typically for percolation two distinct criteria can be 
adopted to define the critical occupation fraction p c [l2j ■ 
Very often the critical point is defined as the largest value 
of the bond occupation probability p where the largest 
cluster occupies a vanishing fraction of the system in the 
thermodynamic limit, L — > oo. Alternatively, the critical 
point can be defined as the smallest occupation fraction 
for which exists a cluster connecting opposite sides of the 
lattice. In the traditional percolation model, q = 1, these 
two definitions are equivalent For clarity, since we 
halt the occupation of bonds when a spanning cluster is 
formed, we are here examining the behavior of the system 
in terms of the later definition. The scaling relations are 
used to identify the critical condition. Moreover, the crit- 
ical fraction and the exponents resulting from our simu- 
lations for q — 2 are consistent with the values reported 
in Ref. [4| for the Achlioptas model on square lattices. 

In Fig. Q] we show typical clusters built with the gen- 
eralized EP model for q = 1 (traditional percolation 
model), q — 2 (Achlioptas model (H), and q = 10. As 
depicted in Fig. QJi, the case q — 1 results in a standard 
spanning cluster with fractal dimension d c i us = 1.89 
which includes a conducting backbone of significantly 
smaller dimension as a subset, di, ac k = 1.64 [l2j, and a 
few cutting bonds, d cu t = 0.75. In the Achlioptas model, 
q = 2, the spanning cluster incorporates a much larger 
fraction of the system, while the conducting backbone be- 
comes more tenuous with cutting bonds appearing more 
frequently. As q increases, both spanning cluster and 
conducting backbone become larger and more compact, 



occupying most of the lattice, with the number of cut- 
ting bonds being substantially reduced. Surprisingly, the 
plots of the sample clusters shown in Fig. Q] suggest that 
some aspects of the system geometry (e.g., the conduct- 
ing backbone) do not necessarily change monotonically 
as q increases from one to very large values. 

In order to characterize the critical properties of the 
model, we perform finite-size scaling analysis. Figure [2] 
shows the average values of the quantities M r j us , M{, a ck, 
and M cut at criticality against the system size L for q = 1 , 
4, 16 and 64. The scaling behavior with size for the con- 
ductivity a is shown in Fig.|3]for the same values of q. As 
depicted, for any value of q, every measure depends on 
I as a power-law. In Fig. [4] we show that the fractal di- 
mension of the spanning cluster increases monotonically 
from ddus = 1.89 at q = 1 to d c i us — 2.0 at large q values. 
The fractal dimension of the conducting backbone goes 
from dback = 1-64 for q = 1 to db ac k = 2.0 at large val- 
ues of g, but passing through a minimum d^ack = 1-52 at 
q = 2. On its turn, the fractal dimension associated with 
the number (mass) of cutting bonds initially grows from 
dcut — 0.75, reaches a maximum value of d cu t = 1.06 at 
q = 4, and then decreases to values close to zero for suf- 
ficiently large values of q. Similar to the exponent of the 
cutting bonds, the exponent [i associated with the con- 
ductance also displays a maximum around q = 4. Start- 
ing from the standard percolation value of /j, = 0.98 at 
q = 1, it reaches /it = 1.18 at g = 4, and then decreases 
monotonically, approaching fi = in the limit of large q 
values. For completeness, the inset of Fig. |4] also shows 
the behavior of the critical occupation fraction as it ap- 
proaches the limit p c = 1 as q grows. 
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FIG. 2: (color online) Size dependence of the number of cut- 
ting bonds, Mcut, the number of sites in the conducting back- 
bone, M;, ac fc, and the number of sites in the spanning cluster, 
M c i us . All quantities are averaged over at least 2000 realiza- 
tions precisely at the instant in which a percolation cluster 
is formed. For comparison, we present results for traditional 
bond percolation q — 1, as well as for q — 4, 16 and 64. The 
error bars are smaller than the symbols. In all cases, the mea- 
sured quantities show a typical power-law dependence with 
the system size L (solid lines). The corresponding values of 
the slopes of these curves are shown in Fig. [4] as a function of 



When the occupied bonds are selected from a large 
set, q — > oo, the behavior of the spanning clusters is 
equivalent to the one observed in Ref. [llj]. Since the 
selection rule favors the occupation of bonds connecting 
small clusters, the growth of the spanning cluster is sup- 
pressed, leading to a significant increase in the critical 
occupation fraction p Cl as compared to the standard per- 
colation value p c = 0.5. For very large values of q, the 
spanning cluster is rather compact, saturated with inter- 
nal bonds. All critical exponents obtained are consistent 
with homogeneous networks, namely, d c i us = dback = 2, 
with the critical fraction converging to its maximum pos- 
sible value, namely, p c = 1. This indicates that, for very 
large values of q, the transition is similar to percolation in 
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FIG. 3: (color online) Variation with size of the global con- 
ductance <r, for q = 1, 4, 16 and 64. As for the geometrical 
properties shown in Fig. [21 each point represents an average 
over at least 2000 realizations precisely at the critical per- 
colation point. In all cases, the error bars are smaller than 
the symbols. The solid lines are the least-squares fits to the 
simulation data of power-laws, a ~ Z/~ M . The corresponding 
exponents fi are shown in Fig. [4] together with the results for 
other values of q. 



one dimension. Surprisingly, the behavior at intermedi- 
ate values of q is the least favorable for conduction, that 
is, the conductance decreases faster with system size. At 
this condition, one also observes that the fractal dimen- 
sion of the backbone dback is minimal and the exponent 
controlling the scaling of the cutting bonds is maximal. 
As we discuss next, the reason for this behavior is closely 
related to the product rule selection criteria as well as 
the degree of heterogeneity of the cluster size distribu- 
tion prior to the critical point. 

During the EP process, bonds selected from the g-set 
can be either merging bonds, i.e., bonds that connect 
sites belonging to two distinct clusters, or internal bonds, 
i.e., bonds connecting two sites inside the same cluster. 
Note that only the inclusion of internal bonds can form 
loops in the system. For any value of q > 1, however, the 
inclusion of internal bonds in the larger clusters is hin- 
dered. Larger clusters have more internal bonds than 
smaller clusters, so it is more likely that the internal 
bonds included in the selection set already belong to the 
largest clusters in the system. However, since these bonds 
have the largest possible weight, their occupation is less 
probable. Now we argue that the fractal dimension of 
the backbone is dominated by the proportion of loops. 
For very large values of q, all clusters have more or less 
the same size and, at p c , two nearly compact clusters are 
joined to form the spanning cluster, which means that a 
considerable fraction of the resulting backbone consists 
of loops. Therefore, the fractal dimension of the back- 
bone for large q is close to that of the spanning cluster 
itself, both approaching the topological dimension of the 
lattice in the limit q — > oo. For intermediate values of q, 
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FIG. 4: (color online) The behavior of characteristic scaling 
exponents obtained for the quantities investigated in Figs, m 
and[3] As shown, the fractal dimension of the spanning cluster 
ddus grows with q, approaching d c i us — 2, already at moder- 
ate values of q. In this way, the increase of the parameter q 
leads to a percolation transition that is more abrupt. Other 
exponents display an interesting non-monotonic dependence 
on q. For instance, the fractal dimension of the conducting 
backbone dback is minimal at q = 2, while the exponents for 
the conductance [i and number of cutting bonds d cut show a 
maximum around q = 4. At large values of q, the conducting 
backbone becomes denser, dback = 2, with a conductance that 
is practically independent on the system size, p ~ 0. These 
results point to a first-order transition where percolation clus- 
ters are compact and occupy a large fraction of the lattice. 
Note also that the critical connectivity increases with q ap- 
proaching p c — 1 in the limit of large values of q, as shown in 
the inset. 



however, the cluster size distribution is not that strongly 
monodisperse and, due to the selection rule, the connec- 
tion of clusters of intermediate sizes dominates over the 
closing of loops in the larger clusters. At p c the spanning 
cluster is typically created by joining two of the larger 
clusters. Thus the backbone will have a particularly low 
proportion of loops and therefore its fractal dimension 
will be low. This explains the observed minimum. 

In summary, we have studied the transport properties 
of percolation networks built from a generalization of the 
Achlioptas model 0, where a bond is occupied from a 
set of q > 1 randomly selected bonds, according to the 
product rule of cluster sizes that they would potentially 
connect. Our results indicate that, for any value of q, 
the mass of the conducting backbone, the total mass of 
cutting bonds, and the global conductance of networks 
at the critical point, all scale as power- laws with system 
size. Moreover, the scaling exponents converge to the ex- 
ponents of a fully occupied square lattice, in the limit of 
very large values of q. This behavior is consistent with 
a first-order phase transition. Interestingly, we observe 
that systems with intermediate values of q display the 
worst conductive performance. We argue that conduc- 
tion becomes difficult for these values of q because the 
formation of loops in the spanning cluster is prevented, 
resulting in smaller conducting backbones. 
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